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Abstract 

A continuum (Mullins-type) model is formulated for the isotropic evolution of a 
solid surface on which the mass transport occurs by oscillatory surface diffusion. The 
time-space oscillations of diffusivity are assumed to be induced by external source(s); 
a number of possible sources, different in their physical nature has been suggested in 
[1]. The particular driving mechanism assumed here is low-energy, pulsed electron or 
laser beam(s) incident on a surface. The simplified version of the model is applied to 
study the relaxation of a periodic surface corrugation. Several extensions of the model 
are proposed. 

Keywords: Surface diffusion, Surface structure, morphology, roughness, and topog- 
raphy, Semiconducting surfaces, Curved surfaces, Models of non-linear phenomena. 

Considered in this Letter is a model 1+1 case corresponding to a two-dimensional crystal 
with a one- dimensional surface. The evolution of the crystal surface given by z = h(x, t) is 
described by the following phenomenological partial differential equation [2, 3]: 

;,/ V = -<>./,. (1) 
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where V is the normal velocity of the curve representing the surface, Q = const, the atomic 
volume, j the surface diffusion flux, and s the arc length along the curve. Subscripts denote 
differentiation in all equations in the Letter. 
The surface flux j is given by 

Dv 

where D is the surface diffusivity of adatoms, v = const, the surface density of adatoms, k 
the Boltzmann constant, T the temperature, and 

a = Qjk (3) 
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the surface chemical potential. In (3), 7 is the isotropic surface energy density and 



the surface mean curvature. Note that more complicated expressions for [i must be used in 
the cases of weakly and strongly anisotropic surface energy density [4] . 

1 Oscillatory driving of surface diffusion by pulsed elec- 
tron or laser beams 

It was recently pointed out in [1] that "oscillatory driving appears as an alternative route for 
pattern formation, with two basic advantages: First, patterning and growth are separated, 
so that morphology is not a function of the growth process. Second, it offers better control 
of the structure. An in situ and real-time control of the pattern becomes possible, opening 
a wide range of new applications." On the side of fundamental physics, oscillatory driving 
of solid surfaces falls in the scope of intensively studied nonlinear phenomena such as the 
inverted pendulum problem, the Faraday instability of a fluid surface, and pattern formation 
in granular media (see references in [1]). 

Work [1] employs a mesoscopic step-flow model for the surface to calculate the induced 
mass flux along vicinal surfaces, which accounts for surface dynamics and stability. Slope 
selection, surface metastability, and driving frequency-dependent surface stability are found. 
No particular physical driving mechanism was assumed in [1], and only the surface diffusivity 
oscillating in time but uniform in space was considered. 

Let pulsed, low-energy photon or, more often, electron beam is incident on a solid surface 
('low-energy' means that the incident energy is insufficient to cause the removal of atoms 
from the surface, through elastic collision). The collision energy generates an excited surface 
state, often free carriers (electrons and holes), which recombine with the local release of 
heat. Two types of enhanced surface diffusion are usually distinguished ([5], pp. 250-254). 
One mechanism is 'local heating' or 'hot-spot'; the key step here is the conversion of the 
recombination energy into local vibration energy of atoms. Some of the recombination energy 
and some thermal energy are used to achieve a passage over the potential barrier to motion. 
Another mechanism is 'local excitation' and it makes use of an energy surface corresponding 
to an excited electronic state. In reality, at least when laser beams (photons) are incident on 
a semiconductor surface both mechanisms are accompanied by 'direct' enhancement, that is 
by conversion of incident energy into lattice vibrational energy through collisional momentum 
transfer. 

In most cases, one recombination event can produce at most one diffusion jump. This 
means that rather large number of recombinations is needed to cause significant matter 
transport. For the atomic displacements through collisional momentum transfer, the energy 
of electrons should be higher than about 100 keV; for displacement through electronic ex- 
citation, the electron energy could be as small as 5 eV (9 eV for GaAs; see Table 2.1, p. 
60 in [5]). Publications [6]-[ll] discuss surface effects caused by electron and laser beams 
(including the diffusion). 



-h 



XX 



(4) 



K = 




2 



On a macroscopic level of description used here, pulsed radiation gives rise to a quasista- 
tionary state in which the temperature at the targeted spot fluctuates about the mean value 
T with a frequency equal to the source pulse repetition frequency u, and an amplitude T(x), 
where x is coordinate across the spot [12]. The value of T and the functional form of T(x) 
are determined by the impulsive power density / and the mean power density I = utj of 
the radiation (tj is the pulse duration), the absorbivity of the surface at the radiation wave- 
length, and the thermophysical characteristics of the material. Thus the surface diffusion 
can be, at least in principle, controlled through temperature oscillations, and I consider this 
mechanism in the following sections. 



2 Derivation of the surface evolution equation 

Let the surface diffusivity is given by the familiar Arrhenius expression 



D(T) = A) exp 
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where D = const, and E d = const, are the pre-factor and activation energy, respectively. 

Substituting the equations (2) - (5) into the equation (1), and using d/ds — (1 + 
h 2 x )~ l l 2 d/dx I obtain the following 4th-order surface evolution equation: 
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where B = Q 2 uD '~f. 

The equation (6) admits trivial solution (base state) 

h = h = const., T = T Q = const. > 0, 



(7) 



corresponding to the unperturbed, horizontal surface at elevated mean temperature. With- 
out loss of generality ho — can be chosen. Consider small perturbations of this state: 

h = ~h{x,t),T = T + f{x,t), where \h(x,t)\, \f(x,t)\ = ^ < 1. (8) 

To 

One can think of a number of cases for the perturbation of the mean temperature field, each 
of which may result in very different nonlinear phenomena at the surface. For instance: 
several, synchronized radiation sources are arranged in a periodic array along the surface. 
Then, the spatial shape of the temperature perturbation can roughly be taken periodic with 
the period equal to the separation between the beams. Or, each source has its own pulse 
repetition frequency; then, the resulting situation amounts to a multi-frequency driving of a 
surface (also see Section 4). 

Expanding exp (^f^) / (kT) in (6) in powers of T(x, t) and keeping only the zeroth-order 
and linear term, I find: 



h t = Ba 2 f x {\ + h 2 x ) 
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where 

Next, the equation (9) is expanded in powers of h x up to second order; considered are only 
such perturbations of the temperature field that |!r x | < \h x \ <C 1 (so that terms proportional 
to Th x ,T x h x and higher order are ignored). This procedure yields the following equation: 

3 



h t = B [-(oi + a 2 T)h xxxx - a 2 T x h xxx + 3(ai + a 2 T)h xx + 

10 ai h xx h xxx h x + ai (2/w* - 27/4) &*] + C(^). (11) 

The equation (11) is the starting point for studies of effects of temperature variations on the 
surface diffusion-produced morphologies of solid surfaces. More terms of the expansion must 
be kept in this equation if the condition of small slope of the temperature perturbation is 
relaxed. The key simplifying assumption in the derivation is the independence of the surface 
energy density on temperature, and its isotropy. 



3 Case of spatially uniform temperature field 

Let the temperature perturbation is the time-periodic function, that is T(x, t) = T(t) 
Tcosut, <T<1. The equation (11) becomes 



h t = -B 



(ai + a 2 T cos ut) (h xxxx - 3/4) ~ ^Oaih xx h xxx h x - a x (2h xxxx - 27/4) K 

(12) 

Notice that (12), as well as the more general equation (11) is nonlinear even in the zeroth 
order of small-slope expansion. Simple linear approximation is now considered (all spatial 
derivatives <Cl). 

In the linear approximation, (12) reduces to 

h t = -B(l + Acost)h xxxx , B = ^->0, A=—T. (13) 

Equation (13) is in the nondimensional form; 1/ou is used for unit of time and the diameter of 
the irradiated spot, L, for unit of length. Notice that A is positive if E d > kT and negative 
otherwise. 

Solution of the equation (13) implies that the initial surface corrugation with the shape 
of cos arc (where a is arbitrary wave number) relaxes to the base state h = in oscillatory 
fashion as 

exp (-Ba 4 t) exp (-£Ma 4 sin t) , (14) 

while preserving the initial shape. Such relaxation mode should be contrasted to monotonous 
exponential decay as 

exp (-BaH) (15) 

in the absence of oscillatory driving, that is at constant T = T Q (here, for comparison with 
the case of oscillatory driving, the time and length units are chosen same). At typical values 
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L = 100 nm, Q = 2 x 10" 23 cm 3 , 7 = 10 3 erg/cm 2 , v = 10 15 cm" 2 , to = 5 x 10 4 s _1 , D = 10 12 
cm 2 /s, E d = 0.5 eV, T = 300 K, T = 10~ 2 , one obtains B = Ax 10 4 , A = 0.2. Thus the 
corrugation relaxes fast unless it has long wave length; with the diameter of the irradiated 
spot or pulse repetition frequency increased, shorter wave length perturbations live longer. 

4 Discussion and conclusions 

This Letter suggests the simple, 1+1 continuum model of the driven evolution of the mor- 
phology of a crystal surface irradiated by the pulsed, low-energy electron or laser beam(s). 
The solution to the evolution equation has been obtained analytically in the particular case 
of spatially uniform field of temperature oscillation on a surface, and the small-slope, linear 
approximation. 

A number of future directions for the research can be identified. 

1. Study of the evolution of the surface morphology within the context of nonlinear, small- 
slope equations (12) and (11), as well as of the equation (9) (which is applicable for 
large slopes); the latter two equations require an assumption on the form of spatial 
dependence of the temperature perturbation. These equations are unsolvable analyti- 
cally if no additional simplifying assumptions are made. Therefore, one has to resort 
to numerical methods. It is expected that the nonlinearity will support the oscillations 
and slow down the rate of relaxation of the corrugation with the emergence of quasi- 
stationary pattern. Its structure and dependence on driving frequency and amplitude 
is of prime interest. 

2. Derivation, analysis and solution of the surface evolution equation under assumption 
that temperature oscillations exert influence on the anisotropic surface energy density 
only. Note that during crystal growth with anisotropic surface energy (no oscillatory 
driving) and in the longwave limit a stationary solutions in the form of a periodic hill 
and valley structures exist [13]; the oscillatory driving is expected to produce a different 
quasistationary pattern. The model that is most relevant to practice would include 
the temperature dependence of both the diffusivity and the surface energy density. 

3. The drawback of the model presented here is the ad hoc inclusion of the temperature 
perturbation. To make the model more realistic, one should formulate and solve the 
problem for the thermal diffusion on the irradiated surface [12, 14]. The solution of 
that problem, T(x,t) must then be used for the temperature perturbation. 

4. Modeling of the crystal growth. Assume for a moment that the oscillatory driving 
of a crystal surface through the pulsed radiation mechanism discussed in Section 1 
is absent. Then the impinging oscillatory flux of material particles on a surface will 
not result in new instabilities but may change the thresholds and resulting patterns 
arising from known instabilities during crystal growth with non-oscillatory flux [1]. 
However, the two oscillatory mechanisms acting simultaneously may create completely 
new instabilities. This situation bears resemblance to two-frequency driving of a fluid 
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surface, which is an actively studied subject [15] (of course, physical media, conditions, 
mechanisms and governing equations are completely different). 

5. Extension of the above models to second spatial dimension. 

I'm currently working on the first two directions. The long-term impact of the mathemat- 
ical modeling of externally driven surfaces to the technologically important field of studies 
of surface nanoscale patterns and their applications is expected. 
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